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ABSTRACT 

We present Rc-band surface photometry for 170 of the 203 galaxies in GHASP, Gassendi 
H-Alpha survey of SPirals, a sample of late-type galaxies for which high-resolution Fabry- 
Perot H maps have previously been obtained. Our data set is constructed by new Rc-band 
observations taken at the Observatoire de Haute-Provence (OHP), supplemented with Sloan 
Digital Sky Survey (SDSS) archival data, obtained with the purpose of deriving homogeneous 
photometric profiles and parameters. Our results include Rc-band surface brightness profiles 
for 170 galaxies and ugriz profiles for 108 of these objects. We catalogue several param¬ 
eters of general interest for further reference, such as total magnitude, effective radius and 
isophotal parameters - magnitude, position angle, ellipticity and inclination. We also perform 
a structural decomposition of the surface brightness profiles using a multi-component method 
in order to separate disks from bulges and bars, and to observe the main scaling relations 
involving luminosities, sizes and maximum velocities. 

We determine the Rc-band Tully Fisher relation using maximum velocities derived solely 
from Ha rotation curves for a sample of 80 galaxies, resulting in a slope of —8.1 ± 0.5, 
zero point of —3.0 ± 1.0 and an estimated intrinsic scatter of 0.28 ± 0.07. We note that, 
different from the TF-relation in the near-infrared derived for the same sample, no change in 
the slope of the relation is seen at the low-mass end (for galaxies with Vmax < 125 km/s). 
We suggest that this different behaviour of the Tully Fisher relation (with the optical relation 
being described by a single power-law while the near-infrared by two) may be caused by 
differences in the stellar mass to light ratio for galaxies with Vmax < 125 km/s. 

Key words: galaxies: photometry, galaxies: structure 


1 INTRODUCTION 

Historically, spiral galaxies have performed a critical role in the 
studies of the dark matter. Observations of the outer flat ro tation 
curves in spiral galaxies (e.g. lRubin. Thonnard. & Fordlll978h have 
focused t he attention to the then overlooked missing mass problem 
(see, e.g. lZwic^ll937h . that stresses the fact that most of what we 
see (light) is just a fraction of what we would like to observe (mass). 
A critical further step, yet to be accomplished, is to understand the 
connection between ordinary and dark matter in the inner regions 


* Based on obseivations performed at Observatoire de Haute Provence, 
France 
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of galaxies to understand whether (and possibly how) light traces 
mass. 

The kinematic decomposition of velocity fields of spiral galax¬ 
ies is th e general method to map their distribution of dark mat¬ 
ter (e .g., vmi^lbadaetalj[l^85 ; van Albada & Sancisill 19861 : iKend 


1 19861 : Kassin. de Jong. & Weineill2006h . However, the stellar mass 

distribution is poorly constrained, and the under-determined stellar 
mass-to-light ratio (M/L) translates into degeneracies, such as the 
disk-halo and the cusp-core problems, that prevent unique decom¬ 
positions. In this context, high resolution, accurate rotation curves, 
such as the observed by the Gassendi H-Alpha su rvey of Spi¬ 
rals ( GHASP), are necessary to alleviate the problem dPutton et al.l 
| 200 ^ . 

Previous works have supported the scenario of cored dark mat- 
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ter profiles le.e. ISpano et al.ir2008h . but studies on the systematic 
errors and larger, homogeneous samples, are still needed to con¬ 
firm these results. This series of papers on the GHASP survey has 
the goal of imposing tighter constraints on the study of dark matter 
distributions in spiral galaxies. In this paper, we build a new surface 
photometry data set for 128 GHASP galaxies in the Rc-band, ob¬ 
served over the years at the Observatoire de Haute-Provence (OHP), 
which provides the basis for the determination of stellar masse 
in forthcoming work. Additionally, we complement this data with 
public Sloan Digital Sky Survey (SDSS) data in order to obtain 
ugriz photometry for 108 GHASP galaxies as well as to increase 
the Rc-band data to 170 galaxies (« 84% of the survey). 

Besides the surface brightness profiles, we also compile a ho¬ 
mogeneous photometric catalogue including several photometric 
quantities of general interest, such as magnitudes, sizes and isopho- 
tal properties. In addition, we perform a multi component light de¬ 
composition in order to separate the light from the disks (our main 
interest to dynamical decomposition) from other components such 
as bulges and bars. Finally, we perform a first set of applications 
to our data set by determining important scaling relations with lu¬ 
minosity, size and velocity of galaxies, and by deriving the Tully 
Fisher relation in the Rc-band. 

This paper is organized as follows. The GHASP sample is 
briefly outlined in Section Following this, the details of the ob¬ 
servations, data reduction and calibration are shown in Section]^ 
In Section|4]we present the methods used for the determination of 
the surface brightness, PA, ellipticity and integrated magnitude pro¬ 
files, and we detail the process of multi component decomposition. 
In Section [S] we test our results against other similar works, and 
we check the internal consistency of our results. Finally, in Sec- 
tion[6]we derive several scaling relations involving luminosity, size 
and rotation velocity using the decomposition results, with special 
emphasis on the Rc-band Tully-Fisher relation. 


2 THE GHASP SAMPLE 


The GHASP sample consists of 203 spiral and irregular galaxies in 
the local universe for which high-resolution Ha maps have been 


observed with Fabry-Perot interferometry (Garrido et al 

2002, 

20031: iGarrido. Marcelin. & AmramI |2004|; Garrido et al 

2005: 

Snano et alj|2008l: lEninat et alj 20081: 

Eninat. Amram & Marcelin 

20081: Eninat et alj 2010l: Torres-Flores et al, 2011). The GHASP 


sample was initia lly designed to be a subsample of the Weste rbork 
survey (WHISP. lvan der Hulst, van Albada, & Sancisilll200lh with 
the goal of providing a local universe reference for kinematics and 
dynamics of disk-like galaxies. 

The GHASP sample was designed to cover a large range of 
morphological types, including ordinary, mixed-type and barred 
galaxies, thus excluding only early-type galaxies because of their 
low Ha content, as illustrated in Figure [T] The photometric sam¬ 
ple presented here is built with data coming from two sources. 
Photometric Rc-band data was obtained by the GHASP collabora¬ 
tion at the OHP over the last decade for 128 galaxies. To enlarge 
the sample, we also take advantage of the public dataset from th e 
seventh Data Release (DRV) of the SDSS ( fstoughton et alj|2002h . 
which provides imaging and calibration in five pass bands (ugriz) 
for 108 of our galaxies. By the combination of both data sets, we are 
able to obtain photometry of a total of 170 GHASP galaxies, which 
are listed in table [T] with details about the observation. However, 
we note that, on average, the data observed at OHP in the Rc-band 
goes about half magnitude deeper than the SDSS data. 
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Figure 1. Frequency of galaxy morphological types of the GHASP sample 
according to the Hyperleda classification, separated in sections considering 
only ordinary spirals (59 galaxies), mixed types (52), barred (70), irregulars 
(19) and lenticulars (2). 


Table 2. List of observational runs in which the Rc photometry of the GHASP 
galaxies was obtained at OHP. 


Run 

Dates 

Number of 

Galaxies observed 

(1) 

2002 Mar 7th - Mar 13rd 

29 

(2) 

2002 Oct 28th - Nov 10th 

38 

(3) 

2003 Mar 8th - Mar 9th 

17 

(4) 

2003 Mar 29th - Apr 6th 

22 

(5) 

2003 Sep 22nd - Sep 28th 

25 

(6) 

2003 Oct 21st - Oct 25th 

6 

(V) 

2008 Jun 2nd - Jun 4th 

11 

(8) 

2009 Oct 23rd 

1 

(9) 

2010 Mar 19th-Mar 21st 

16 


3 DATA REDUCTION 

3.1 Data from the OHP observatory 

Broadband imaging for 128 galaxies in the Rc-band was obtained 
with the 1.2m telescope at the Observatoire de Haute-Provence 
(OHP), France, in several observation runs as presented in Table 
12] The images have a field of view of 11.7’ x 11.7’, taken with a 
single CCD with 1024 x 1024 pixels, resulting in a pixel size of 
0.685 arcsec“^. 

Basic data reduction was performed with IRAlQ tasks, includ¬ 
ing flat-field, bias subtraction and cosmic-ray cleaning. Images of 
the same galaxy are then aligned and combined for the cases with 
roughly the same smallest seeing FWHM, estimated from isolated 
field stars. Photometric stability check and zero point calibration 
was obtai ned by the obse rvation of several standard stars from the 
catalog of lLandoi^ ( Il992h in different times during the nights, con¬ 
sidering the mean airmass correction coefficient of 0.145 for the 
Rc-band dChevalier & Ilovaiskvlll99ih . and no colour term. 

The determination of the sky level is the greatest sou rce of un¬ 
certai nty for surface brightness profiles and magnitudes dCourteaui 
Il996h . For this purpose, we adopt the method of estimating the 
background by selecting sky “boxes” on the images, which are vi¬ 
sually selected areas where the galaxy and stellar light contribution 
is minimal, and use those regions to calculate a smooth surface us- 


^ IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in As¬ 
tronomy (AURA) under cooperative agreement with the National Science 
Foundation. 
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Table 1. Photometric sample of the GHASP survey. The printed version contains only an abridged version of the table, with the remaining material available 
online as supplementary material. (1) Galaxy name. (2-3) Right ascension and declination of the galaxies according to the NED database (4-5) Morphological 
classification accor ding to the Hubble type and to the de Vaucouleurs numerical type from the Hyperleda database. (6) Distance to the galaxies according to 
lErrinat et al.l ilOOSi) . (7) Shows if the galaxy is in the SDSS. (8-10) OHP Rc-band observation log, including the runs, total exposure time and the seeing. 


Galaxy 

(1) 

a 

(J2000) 

(2) 

5 

(J2000) 

(3) 

Morphology 

(4) 

Morphology 

t 

(5) 

Distance 

(Mpc) 

(6) 

SDSS 

(7) 

OHP observation log. 

Runs Exptime FWHM 

(s) (arcsec) 

(8) (9) (10) 

UGC 89 

00h09m53.4s 

+25d55m26s 

SBa 

1.2 ±0.6 

64.2 

no 

2,5 

600 

2.7 

UGC94 

00hl0m25.9s 

+25d49m55s 

S(r)ab 

2.4 ±0.6 

64.2 

no 

5,6 

3300 

2.3 

IC 476 

07h47ml6.3s 

+26d57m03s 

SABb 

4.2 ±2.6 

63.9 

yes 

— 

— 

— 

UGC 508 

00h49m47.8s 

+32dl6m40s 

SBab 

1.5 ±0.9 

63.8 

no 

5 

3600 

2.3 

UGC 528 

00h52m04.3s 

+47d33m02s 

SABb 

2.9± 1.1 

12.1 

no 

2 

1500 

1.9 

NGC 542 

01h26m30.9s 

+34d40m31s 

Sb pec 

2.8 ±3.9 

63.7 

yes 

— 

— 

— 

UGC 763 

01hl2m55.7s 

+00d58m54s 

SABm 

8.6 ± 1.0 

12.7 

yes 

2,5 

600 

3.1 

UGC 1013 

01h26m21.8s 

+34d42mlls 

SB(r)b pec 

3.1 ±0.2 

70.8 

yes 

2,5 

5100 

2.5 

UGC 1117 

01h33m50.9s 

+30d39m37s 

Sc 

6.0 ±0.4 

0.9 

no 

5 

4500 

2.7 

UGC 1249 

01h47m29.9s 

+27d20m00s 

SBm pec 

8.8 ±0.6 

7.2 

no 

2 

1800 

2.1 

UGC 11951 

22hl2m30.1s 

+45dl9m42s 

SBa 

1.1 ±0.8 

17.4 

no 

2 

2100 

2.6 

UGC 12060 

22h30m34.0s 

+33d49mlls 

IB 

9.9 ±0.5 

15.7 

no 

6 

6000 

2.5 

UGC 12082 

22h34ml0.8s 

+32d51m38s 

SABm 

8.7 ±0.8 

10.1 

no 

6 

3000 

3.7 

UGC 12101 

22h36m03.4s 

+33d56m53s 

Scd 

6.6 ±0.9 

15.1 

yes 

2 

1800 

1.9 

UGC 12212 

22h50m30.3s 

+29d08ml8s 

Sm 

8.7 ±0.5 

15.5 

yes 

2 

1800 

2.2 

UGC 12276 

22h58m32.5s 

+35d48m09s 

SB(r)a 

1.1 ±0.5 

77.8 

no 

2 

2700 

2.0 

UGC 12276c 

22h58m32.5s 

+35d48m09s 

S? 

5.1 ±5.0 

77.8 

no 

2 

2700 

2.0 

UGC 12343 

23h04m56.7s 

+12dl9m22s 

SBbc 

4.4 ±0.9 

26.9 

no 

2,5 

1500 

2.7 

UGC 12632 

23h29m58.7s 

+40d59m25s 

SABm 

8.7 ±0.5 

8.0 

no 

5 

4500 

2.4 

UGC 12754 

23h43m54.4s 

+26d04m32s 

SBc 

6.0 ±0.4 

8.9 

no 

2 

1200 

2.3 


ing the IRAF package IMSURFIT with polynomials of order 2, which 
is subtracted from the original images. This process resulted in an 
homogeneous background for which the typical residual standard 
deviation is in the range 0.5-1% of the sky level. 

We have modeled the Point Spread Function (PSF) of our im¬ 
ages using the IRAF PSFMEASURE task. We selected bright, un¬ 
saturated stars across the fields using the task DAOEIND, a nd then 
mode led their light profiles using a circular Moffat function ( iMoffall 
1 19691) . given by 


method for Rcimages, and also because a few authors have pointed 
out errors in the sky deter minations on images with bright galaxies 
in early SPSS rele ases (e.g. lSemardi et ai]|2007l : lLauer et al.ll2007l : 
iLisker et alj2007l) . 


4 DATA ANALYSIS 
4.1 Surface photometry 


PSF(r) 



( 1 ) 


where the radial scale length a and the slope (3 are free parame¬ 
ters, which can be relate d to the seeing by th e relation FWFIM= 
2a\/2'^/P — 1 (see also iTruiillo et al.1 l200ll) . This method has 
proved to be suitable in our case due to the presence of extended 
wings in the PSFs. Overall, the typical seeing of our observations 
is FWHM«3 arcsec, with the parameters a, /3 and FWFIM having 
mean statistical uncertainties of 1.7%, 9% and 3% respectively. 


3.2 Data from SDSS 

To increase the number of gal axies in our photometric sample in the 
Rc-band, we use SDSS DR7 dAbazaiian et al.ll200^ archival data 
for 108 GHASP galaxies we found in the database and transformed 
SDSS ugriz data into RcWith a multi-band scaling relation (more 
details in section IdTt . Among these 108 galaxies, 66 have also been 
observed in the OFIP, thus 42 new ones are added to the final photo¬ 
metric sample. Calibrations and the PSF of the images are obtained 
directly from the data products of the survey. We performed a new 
sky determination for each image for consistency with the adopted 


We study the photometric properties of th e sample using the tradi¬ 
tional method of elliptical isophote fitting ( lKentil984l : ljedrzeiewskil 
Il987h . Surface brightness (SB) profiles of the galaxies were ob¬ 
tained using the IRAF task ELLIPSE, which provides a number of 
parameters that describe the light of the galaxy as a function of the 
semi-major axis (which we simply refer to as the radius, r), includ¬ 
ing the ellipticity (e), position angle (PA) and the curve of growth, 
which quantifies the total apparent magnitude inside each isophote. 

Masks for foreground and background objects were produced 
interactively in two steps. Firstly, most objects in the images were 
detec ted and masked out with S EXTRACTOR dBertin and ArnoutsI 
1 19961) . Other important sources not detected by the program, such 
as saturated stars and stars/galaxies superposed to the galaxies 
of interest were then masked during ELLIPSE runs. Finally, we 
checked the results inspecting the residual image produced by sub¬ 
tracting an interpolated model of the galaxy produced with the 
task BMODEL. This process was carried out several times for each 
galaxy until no bright sources were observed in the resulting sub¬ 
tracted images except for spiral arms and/or bars of the galaxy not 
masked on purpose. 

The centre of each galaxy was defined in a first iteration of 
ELLIPSE and later was set fixed for all other iterations. The posi¬ 
tion angle and ellipticity of the isophotes were usually set free to 
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vary as a function of the rad ial distance, as usually done for late- 


type g 

alaxy photometry (e.g. 

Balcells et alJl2003l;lMac Arthur et al.l 

|2()Q3l; 

McDonald. Courteau. & Tu11y|| 2009|). but in a few cases we 


were forced to fix the geometric parameters in part or in the whole 
galaxy in order to obtain convergence for the photometry. For the 
SDSS data, we adopt this method in the r-band images, due to its rel¬ 
atively high signal-to-noise ratio in the ugriz system, but we fix the 
position angles and ellipticities accordingly to the r-band parame¬ 
ters in the other pass bands in order to obtain consistent colours. 

Uncertainties for the SB profiles include the isophote deter¬ 
mination error given by ELLIPSE, the photon counting statistics of 
the detector, and the sky level subtraction uncertainty, all added in 
quadrature. All profiles are corrected for t he Galactic foregroun d 
extinction using the dust reddening maps of ISchlegel et al.l ( Il998ll . 
assuming a dust model with constant selective extinction of 3.1, 
and rela tive extinction for the different pass bands according to ta¬ 
ble 6 of ISchlegel et al.l (Il998h . However, we do not attempt for a 
correction of the SB profiles for the more uncertain problem of the 
galaxies’ internal extinction. 

Finally, to obtain Rc-band SB profiles from SDSS data, we use 
a sligh tly modified version of fhe relafion derived by Ijesfer et al.l 
ll20o‘5h . given by 

Mfi(f) = 0.42^g(r) — 0.38^r(r) -f 0.96/ii(r) — 0.16, (2) 

where fi{r) represents the surface brightness profile al radius r in 
the pass band indicated in the subscripts. The equation above was 
originally derived for stellar photometry, so we have tested its ac¬ 
curacy in surface photometry by the comparison of OHP SB pro¬ 
files, obtained directly in the Rc-band, with profiles derived from 
the SDSS ugriz bands using a sample of 54 galaxies for which the 
geometric parameters of both data sets are similar. The results are 
presented in Figure which shows the difference of the profiles 
as a funclion of the OHP surface brightness profiles relative to the 
sky level, which varies in the OHP observations, from one galaxy 
observation to another. The red line shows the running RMS differ¬ 
ence between the profiles, indicating that the error in the process of 
transforming between the photometric systems is of ~ 0.08 mag 
arcsec“^ in the regions brighter than the sky level, ~ 0.15 mag 
arcsec“^ for the regions down to two magnitudes fainter than the 
sky, and ~ 0.4 mag arcsec“^ for the regions 5 magnitudes fainter 
than the sky level. 

We present a sample of SB profiles for a variety of morpho¬ 
logical types in Figure including also the ellipticity and posi¬ 
tion angle variations. All surface brightness profiles are available 
in electronic format. In the next section, we detail other catalogued 
Rc-band photometric properties derived from SB profiles in this 
section. 



Figure 2. Difference between Resurface brightness profiles for OHP and 
SDSS data for 54 galaxies in common to both samples as a function of 
the OHP SB profile relative to the sky level. Profiles for the SDSS galaxies 
are obtained using equation G). Each blue line represents the difference in 
the surface brightness profiles of a single galaxy profile comparison. The 
dashed red lines indicates the running RMS difference between the profiles. 


(PAiso), ellipticity (eiso) and apparent magnitude (mR,iso) directly 
from the SB profiles, with uncertainties estimated by Monte Carlo 
simulations of perturbations of the profile according to their uncer¬ 
tainties. Also, the inclination of the gala xies is estimated at a given 
isophotal level as dTullv & Fisheilll98'ill 


COS iiso 


(1 £iso) Qo 


(3) 


where go = 0.2 is the intrinsic flattening of edge-on disks (e.g. 
iHavnes & GiovaneliHl 19841 : ICourteajll99m . 

We also measured the total (asymptotic) apparent magnitudes 
of the galaxies (ma^totai), which were calculated by the extrapo¬ 
lation of the curv e of growth of the SB profiles using a derivative 
method similar to lCairos et al.l (1200 ih . However, this method failed 
in cases of galaxies for which the curve of growth did not converge. 
In these cases, we used the last isophote total magnitude to estimate 
a lower limit to the total magnitude. Also using the curve of growth, 
we measured the effective radius of the galaxies (rso), which is de¬ 
fined as the radius containing 50% of the total light of the galaxy. 
In the cases where we have not obtained a safe total magnitude, we 
then estimated the lower limits of the effective radius. Uncertainties 
in these parameters are also based on Monte Carlo simulations. 

Table[3]presents a sample of the results for the isophotal level 
of 22.5 mag arcsec'^. The complete catalogue, and the catalogue 
for the isophotal level of 23.5 mag arcsec'^ are provided in the sup¬ 
plementary material. 


4.2 Integrated and isophotal photometry in the Rc-hand 

For the Rc-band SB profiles derived in this work, we obtained a 
number of properties of the galaxies which are of general interest 
by fixing a reference isophofal level. In the case of the Rc-band, the 
isophotal level of 23.5 mag arcsec'^ is usually used as reference, 
because it corresponds to an aperture similar to the B-band isophote 
of 25 mag arcsec'^. However, this level was reached for only 72% 
of our surface brightness profiles. Therefore, in order to provide a 
more complete catalogue for our sample, we also use the isophotal 
level of 22.5 mag arcsec'^ to provide parameters for 98% of the 
sample. 

We measured the isophotal radius (riso), position angle 


4.3 Multicomponent Decomposition 

In our forthcoming work (Pineda et al. in preparation), we plan to 
study the kinematic properties of a subsample of GHASP galaxies 
with specific photometric properties depending on the relative im¬ 
portance of the disks in comparison with bulges and bars. In order 
to separate the SB brightness profiles into different structural com¬ 
ponents, we proceed to a multicomponent decomposition of the SB 
profiles. 

For that purpose, we use a parametric profile fitting method 
which includes as many components as necessary to separate the 
photometric components - disks, bulges, bars, spiral arms, lens and 
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Figure 3. Examples of surface brightness profiles in the Rc-band, position angle and ellipticity profiles for six galaxies for which data are available both for 
the OHP and for the SDSS data sets. The surface brightness profiles are presented in different colours for each band according to the upper right labels, and 
uncertainties are not presented for the sake of clarity. The middle and lower panels of each galaxy show the position angle and ellipticity profiles respectively, 
with gray representing the r- band results and blac k representing the Rc-band geometry. All profiles are con'ected for the Galactic foreground extinction 
according to the dust maps of ISchlegel et al.l il998h . 


nuclear sources. Ideally, one could use a complete 2D fitting to bet¬ 
ter describe the non axisymmetric components, like the bar, but the 
reliability of the ID method to recover the struct ural parameters is 
comp arable to 2D, at least for the disk component dMac Arthur et al.l 
l2003ll . and also allows the estimation of the integrated properties 
with good accuracy. 

We have performed the decomposition in all Rc-band pro¬ 
files in the GHASP OHP sample. We developed a Python routine 
which performs a weighted chi-square minimization between the 
dat a and a model us ing the Levenberg-Marquardt algorithm (see, 
eg. lPress et alJll992ll . using a PSF convolution of the models with 
a Moffat function (see section [JTt . The input model is set man¬ 
ually according to the observation of photometric features in the 
SB profile and the images of the galaxies. Also, the observation of 
the varying ellipticity and position angle as a function of the ra¬ 
dius usually hinted for the different structural sub components of a 
galaxy. Figure|4]presents examples of structural decomposition of a 
sample of twelve galaxies, illustrating the variety of morphologies 


we have in our sample and the varied decomposition components 
we included. The structural parameters for all galaxies are listed 
in Table |4] In the next section, we give some details on the used 
parametrizations for the different components. 


4.3.1 Disks 

Since early works of lPatters^ ( Il940h and ide Vaucouleurll ( Il958h . 
the intensity profile of disks have been mostly described by a simple 
exponential law, 

Id{r) = , (4) 

where Jo is the central (r = 0) intensity of the disk and h is the 
disk scale length. Usually, we refer to the central intensity in terms 
of surface brightness using the relation /io = —2.5 log Jo- For the 
case of exponential disks, the total apparent magnitude is given by 
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Table 3. Isophotal and integrated photometric parameters in Rcband. The printed version contains only an abridged version of the table, and the remaining 
material is available online. (1) Galaxy name. (2) Data source. (3-7) Position angle, ellipticity, inclination, radius and apparent integrated magnitude at the 
isophote of 22.5 mag arcsec“^. (8) Effective radius. (9) Total apparent magnitude. 


Galaxy 

(1) 

Data 

(2) 

PA 22.5 

(arcsec) 

(3) 

£ 22.5 

(degree) 

(4) 

^22.5 

(degree) 

(5) 

^*22.5 
(arcsec) 
(6) 

m.R,22.5 

(mag) 

(7) 

^50 

(arcsec) 

(8) 

^R,total 

(mag) 

(9) 

UGC 89 

OHP 

-27± 11 

0.22 ±0.03 

40 ±3 

41 ±2 

11.60 ±0.03 

19 ±5 

11.4 ±0.1 

UGC94 

OHP 

-84 ± 2 

0.30 ±0.02 

46 ± 1 

28 ± 1 

13.12 ± 0.03 

16 ± 7 

12.7 ±0.4 

IC476 

SDSS 

-79 ±3 

0.4±0.1 

52 ±9 

16 ±2 

14.7 ±0.1 

> 10.0 

< 14.53 

UGC 508 

OHP 

-61 ±4 

0.11 ±0.05 

28 ±7 

67 ± 1 

11.17 ±0.03 

40 ±1 

10.91 ±0.08 

UGC 528 

OHP 

52 ±2 

0.03 ±0.02 

13 ±3 

65 ±3 

10.16 ±0.01 

21 ± 1 

10.10 ±0.04 

NGC 542 

SDSS 

-34 ± 2 

0.73 ±0.01 

79.1 ± 0.9 

23 ±2 

14.39 ±0.05 

12 ±6 

14.2 ±0.6 

UGC 763 

OHP 

-75 ± 23 

0.12 ±0.02 

29 ±3 

57 ±4 

11.61 ±0.06 

47 ±5 

11.1±0.2 

UGC 763 

SDSS 

-79 ± 35 

0.12 ±0.08 

29± 10 

56 ±6 

11.7 ±0.1 

> 39.0 

< 11.46 

UGC 1013 

OHP 

80 ±2 

0.62 ±0.01 

70.6 ±0.7 

53 ± 1 

12.03 ± 0.01 

27 ±4 

11.7±0.1 

UGC 1013 

SDSS 

80 ±3 

0.64 ±0.01 

72.2 ±0.8 

51 ±5 

12.17 ±0.04 

23 ±3 

11.9±0.1 

UGC 12082 

OHP 

18 ±2 

0.42 ±0.01 

56.1 ±0.9 

19 ±3 

15.0 ±0.2 

49 ±23 

12.8 ±0.5 

UGC 12101 

OHP 

-50 ±2 

0.53 ±0.01 

64.2 ± 0.8 

51 ±2 

12.50 ±0.03 

29 ±2 

12.32 ±0.09 

UGC 12101 

SDSS 

-50 ±3 

0.52 ±0.02 

64 ±1 

49 ±4 

12.66 ±0.06 

28 ±21 

12.3 ±0.9 

UGC 12212 

OHP 

-80 ±5 

0.36 ±0.02 

52 ± 1 

14 ± 1 

15.8 ±0.2 

27 ±26 

14 ± 2 

UGC 12212 

SDSS 

-71 ± 2 

0.3 ±0.01 

46.8 ±0.8 

12 ±5 

15.9 ± 0.6 

> 16.0 

< 14.77 

UGC 12276 

OHP 

-47 ±6 

0.23 ±0.01 

40.8 ±0.9 

30 ±1 

13.00 ±0.03 

19 ±5 

12.6 ±0.2 

UGC 12276c 

OHP 

82 ±2 

0.45 ±0.01 

58.3 ±0.7 

7± 1 

17.0 ±0.2 

> 5.0 

< 16.62 

UGC 12343 

OHP 

38 ±2 

0.30 ±0.03 

47 ±2 

102 ± 1 

10.33 ±0.03 

71 ±9 

10.2 ±0.2 

UGC 12632 

OHP 

22 ±2 

0.53 ±0.05 

64 ±4 

24 ±5 

14.8 ±0.3 

66 ±59 

13 ± 1 

UGC 12754 

OHP 

-13 ±3 

0.29 ±0.01 

46 ±1 

89 ±3 

11.16 ±0.03 

47 ±2 

10.9 ±0.1 


Table 4. Decomposition parameters for the GHASP sample in the Rcband obtained at the OHP observatory. The printed version contains only a sample of 
the table, and the remaining material is available online. (1) Galaxy name. (2-4) Parameters for the Sersic function of bulges, bars and other components. (5) 
Visual classification of the components. (6-10) P arameters for the di sks according to the broken exponential function. (11) Classification of the disks regarding 
to the type of breaks, according to the scheme of lErwin et al.l i2Q08l see text for details). (12) Magnitude of the central point source. 





Sersic 




Disk 




PS 

Galaxy 

Mb 

re 

n 

type 

Md 

hi 

ho 

rb 

a 

Type 

TTTps 


(mag arcsec“^) 

(arcsec) 




(mag arcsec“^) 

(arcsec) 

(arcsec) 

(arcsec) 


(mag) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

UGC 89 

16.7 ± 0.1 

2.3 ± 0.1 

1.0 zb 0.1 

bulge 

21.5 zb 1.0 

35 zb 25 

— 

— 

— 

Type I 

— 


20.3 ± 0.2 

19 ± 1 

0.2 zb 0.1 

bar 

— 

— 

— 

— 

— 


— 


23 ± 1 

29 ± 5 

0.11 zb 0.07 

bar 

— 

— 

— 

— 

— 


— 

UGC 94 

18.7 ± 0.3 

0.8 ± 0.3 

1.6 zb 0.3 

bulge 

20.3 zb 0.1 

14.0 zb 0.8 

— 

— 

— 

Type I 

— 


20.79 ± 0.09 

8.3 ± 0.4 

0.14 zb 0.09 

bar 

— 

— 

— 

— 

— 


— 

UGC 508 

18.6 ± 0.2 

3.4 ± 0.1 

2.65 zb 0.04 

bulge 

19.60 zb 0.04 

26.2 zb 0.2 

— 

— 

— 

Type I 

— 


23.69 ± 0.08 

20.2 ± 0.3 

0.05 zb 0.05 

bar 

— 

— 

— 

— 

— 


— 


21.19 ± 0.06 

7.6 ± 0.3 

0.22 zb 0.02 

bar 

— 

— 

— 

— 

— 


— 

UGC 528 

19.0 ± 0.2 

5.4 ± 0.6 

2.7 zb 0.3 

bulge 

18.6 zb 0.2 

19 zb 2 

— 

— 

— 

Type I 

— 


19.8 ± 0.2 

13.1 ± 0.8 

0.1 zb 0.2 

arms 

— 

— 

— 

— 

— 


— 


20.3 ± 0.4 

22 zb 1 

0.15 zb 0.09 

arms 

— 

— 

— 

— 

— 


— 

UGC 763 

28.0 ± 0.3 

311 ± 183 

5.3 zb 0.3 

bulge 

19.99 zb 0.06 

25 zb 2 

— 

— 

— 

Type I 

— 


23 ± 2 

18 ± 9 

0.05 zb 0.06 

bar 

— 

— 

— 

— 

— 


— 

UGC 12276 

21.1 ± 0.1 

5.4 ± 0.2 

2.56 zb 0.01 

bulge 

20.24 zb 0.09 

14.0 zb 0.1 

7.0 zb 0.3 

45.0 zb 0.1 

1.0 zb 0.2 

Type II 

_ 


27.36 ± 0.05 

0.1 ± 0.1 

6.59 zb 0.04 

bar 

— 

— 

— 

— 

— 


— 

UGC 12276c 

22.74 ± 0.09 

0.0 ± 0.1 

2.99 zb 0.08 

bulge 

20.60 zb 0.09 

3.6 zb 0.2 

— 

— 

— 

Type I 

— 

UGC 12343 

20.9 ± 0.4 

12 ± 3 

2.9 zb 0.3 

bulge 

19.2 zb 0.1 

30 zb 3 

— 

— 

— 

Type I 

— 


21.7 ± 0.4 

56 ± 3 

0.2 zb 0.1 

bar 

— 

— 

— 

— 

— 


— 

UGC 12632 

23.1 ± 0.2 

2.9 ± 0.6 

0.05 zb 0.08 

nucleus 

22.1 zb 0.1 

58 zb 4 

— 

— 

— 

Type I 

— 


23.4 ± 0.2 

9.8 ± 0.4 

0.10 zb 0.02 

bulge 

— 

— 

— 

— 

— 


— 

UGC 12754 

20.36 ± 0.07 

11.8 ± 0.4 

0.43 zb 0.04 

bar 

20.6 zb 0.2 

56 zb 7 

— 

— 

— 

Type I 

15.7 zb 0.2 


22.4 ± 0.4 

24 zb 1 

0.05 zb 0.06 

bar 

— 

— 

— 

— 

— 


— 


2-wIoh^^, (5) 

where h/a = 1 — e is the minor-to-major axis ratio of the galaxy. 

However, devi ations from a sim ple exponential disk were 
already noticed by iFreema 3 jl97(]|) . specially in the form of 
truncations or breaks in the inner profiles of galaxies. Later, 
Ivan der Kruit & Searlel ( Il982h also noticed breaks at large radii of 


disks, and more recently deviations at very low sur face bright¬ 
ness h ave been observed, including upward bends (e.g. lErwin et al.l 
l2005h . 

Based on these observations, lErwin etTI ( 1200 8l) proposed a 
reviewed classification of disks in three categories: Type I, sim¬ 
ple disks well described by exponential disks; Type II, disks with 
downward truncations; and Type III, disks wi th upward bends. 
A local census of disk properties performed bv IPohlen & TruiiU^ 


mdisk = —2.5 log 
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Figure 4. Examples of structural decomposition for two galaxies of the GHASP sample in the Rc-band, UGC 89 and UGC 10897. For each galaxy, we present 
seven panels, containing (a) surface brightness profiles and their decomposition components, (b) fitting residuals, (c) position angle profile, (d) ellipticity 
profile, (e) Rc-band image, (f) ELLIPSE model and (g) residual sigma image. 


( l2006lf of late-type galaxies (Sb-Sm) estimated the fraction of Type 
I galaxies to be only 10%, while 60% are classified as Type II 
and 30% are Type III according to this new classification scheme. 
Therefore, an updated profile for the disks is here adopted whenever 
breaks ar e clearly observed , by using broken exponential profiles, 
given by terwin et al.ll2008h 


4.3.2 Other components 

Apart from the disk, several other components are observed, includ¬ 
ing bulges, bars, arms, rings and lenses. We included those compo - 
nents in the decomposition using a Sersic function dSersi 3 ll968ll . 
given by 


Id{r) = SIo exp • {1 + exp [a(r — '*» ', 

(6) 

where Jq is the central intensity of the disk, hi and ho are the inner 
and outer disk scale length respectively, rt, is the break radius, a 
is the sharpness of the disk transition between the inner and outer 
region (where low a means a smooth transition from the inner to 
the outer disk and high a means an abrupt transition), and 5" is a 
scaling factor given by 


5 = [1 + exp(—art)] • (7) 

For the broken disks in our sample, the total luminosities were 
calculated numerically, given that a solution by the integration of 
equation © is beyond the scope of this work. 

In Tabled we include a classification of the disks according to 
our observations. However, it is important to notice that breaks may 
occur at different radial distances, with different physical interpre¬ 
tations: inner breaks {pr ~ 23 mag arcsec“^) may be related to star 
formation, while outer breaks (pr ~ 27 ma g arcsec"^) may indi¬ 
cate a real drop in the stellar mass density dMartm-Navarro et al.l 
l2012h . Therefore, our classifications are restricted to the mean lim¬ 
iting surface brightness of 24.5 mag arcsec“^. 


Ib{r) = le exp 




(8) 


where Ve is the effective scale of the component (for which 50% of 
the light is within re), le is the intensity at the effective radius, 
and n is the Sersic index. The term is not a free parameter, 
but a function of the Sersic index due to the parametrization of 
the function at the effective radius instead of at the centre. In our 
calculations, we adopted the expressions for bn presented in the 
Appendix [Al] of iMacArthur et aO ( l2003l) . In a first order approx¬ 
imation, bn ~ 2n — 0.33, although the error may be considerable 
for n < 0.5. For simplification, we also rescale the effective inten¬ 
sity to surface brightness using the expression p^ = —2.5 log Iq. 
In the case of equation (8ll, t he total magnitude is given by dCiottil 
ll99ll:lMacArthur et al.ll2003l) 


= -2.5 log 


27r/eree*’”nr(2n) b 
6 ^" a 


(9) 


where r(a;) is the complete gamma function of a variable x. The 
Sersic profile is a generalization of other commonly used profile 
functi ons, such as the exponenti al for n = l,a Gaussia n forn = i 
and a ide Vaucouleurd ’s profile dde Vaucouleur dfloTSh for n — 4. 
Moreover, the Sersic index can also be used as an indicator of the 
kind of component that is being observed. For example, in the op- 
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tical and near-infrared wavelengths, iFisher and DrorvI i2010l) have 
shown that n < 2 may indicate a pseudobulge, whereas n > 2 
may indicate a classic bulge for th e spheroidal c omponents. Be¬ 
sides, bars have typically n ~ 0.7 l lGadottill201 ih . In the Table|4l 
where the decomposition results are presented, we include a classi¬ 
fication of the Sersic function components according to the visual 
inspection of images and profiles, such as bulges, bars, lenses and 
spiral arms. 

In 35 galaxies, a nuclear source is also detected, which may 
be related to different physical processes, such as an active nu¬ 
cleus or a stellar concentration. We have tested two approaches for 
parametrizing these components, using either a Sersic component 
or a single delta function with a peak at r = 0. In 12 cases, the 
former approach resulted in a better description of the nucleus, be¬ 
cause they have slightly larger FWHM than that of the modelled 
PSF and/or because of the different shape of the nuclear source 
compared to a star. These components are described as nucleus in 
tables |4] in the column 5. For 23 galaxies, however, the later ap¬ 
proach of using a delta function resulted in a better description of 
the nucleus in these cases. This delta function has only one free 
parameter, the magnitude of the source (rrips), and its profile is of 
a field star which is described as a Moffat function. These point 
source magnitudes are included in the last row of Table|4] 


5 PHOTOMETRIC INTERNAL CONSISTENCY AND 
LITERATURE COMPARISON 

In this section we make a series of tests on our photometric re¬ 
sults to verify their consistency and to compare them with similar 
results in the literature. We have already made an internal consis¬ 
tency check of our SB profiles in Figure|2l where we observed that 
the SB profiles from SDSS data are similar to those observed with 
direct measurements in the Rc-band. Our SB profiles can also be 
compared with those derived for twelve galaxies in common with 
Ide Jong and van der Kruill i ll994) in the Rc-band, as shown in Fig¬ 
ure |5] where black and red lines show the difference between our 
and the literature data for the OHP and SDSS data sets respectively. 
To make a proper comparison, we have fixed the position angle and 
the ellipticity of the galaxies to mimic the method of those authors 
instead of using free position angle and ellipticity as in section|4T| 
We also limited the comparison to regions greater than the seeing 
of our images. 

The most deviant case is UGC 4256, but the internal consis¬ 
tency of our results for two diff erent data sets indicates a possible 
systematic offset in the data of Ide Jong and van der KruitI i ll 9941) 
for this galaxy. The deviation in the outer region of UGC 508 
can be explained by the l imited field of view in the images of 
Ide Jong and van der Kruill i ll9941) which cuts part of the galaxy. 
In this case, it is also possible to see a large variation of the 
isophotes fainter than 21 mag arcsec^ within the observations of 
Ide Jong and van der Kruid Finally, the problematic case of UGC 
10445 for the OHP data can be explained by the relatively short ex¬ 
posure time for this object, including only one image, which affects 
the accuracy of the sky subtraction. Ap art from these remarks, the 
overal l picture is a good agreement with ide Jong and van der KruitI 
i ll9941) . which is in turn in agreement with several other authors in 
the literature. 

In Figurel^ we show the internal consistency of the isophotal 
and integrated photometric parameters derived from the SB profiles 
by the comparison of the results obtained with the OHP and SDSS 
data sets, including position angle, ellipticity, integrated magnitude 



mR,total (mag) - this work 


Figure 7. Comparison of total apparent magnitudes with the litera¬ 
ture. Coloured symbols ind icate the comparison s with data in the Rc- 
band, i ncluding the works oflHeraudeau & SimienHl994),lGil de Paz e t alj 
I 2 OO 3 I) . iKassin. de Jong. & Poggej l2006l) and lAmorin et alj l2007h . Gray 


symbols indicate the comparisons with magnitudes in pass bands sim¬ 
ilar to Rcwith our magnitudes both for the SDSS (circles) and OHP 


(crosses) samples. This later data include magnitudes fromiHickson et al 

198S) 

Tullv et all 

199d). Ijames et all 

.004). ICabrera-Lavers & Garzon 

2004) 

. Tavlor et al 

l2005h. iDovle et al 

2 OO 5 I). iHernandez-Toledo et al 

2007) 

. Noordermeer and van der HulstI 

20071). iThomas et al, l2008l). 

Hemandez-Toledo & Ortepa-F.sbril J2008I). IMatthews & UsonI J2008t) and 


line represents equality bet\veen 

Kiiwattanawong et al.H201lh. The dashed 


measurements. 


and radius at the isophote of 22.5 mag/arcsec^ and also effective 
radius and total magnitude. We also display the mean residual dif¬ 
ference (< A >) and its standard deviation ((t( A)) for each param¬ 
eter, resulting in compatible measurements in both data sets within 
one standard deviation. 

In Figure |7] we compare our total magnitudes with data 
in the literature. There are just a few works in the literature 
for which total magnitudes are measured in the Rc-band, so 
we also include in the figure measurements in similar filters 
in the literature without any additional correction or extrapola¬ 
tion, w h ich are displ a yed as gray symbols, including ^ lly^t_al.l 
lll996h. Ij ames et all 1 2004). ICabrera-L avers & GarzonI 1 2004). 


iTavlor et al. l200 5l)7 Doyle et al. 1 2005). Hemandez-Toledo et al 


1 2007), Thomas et alj~ 20()8l) . i Hemandez-T oledo & Ortega-Esbri 

1 2008), iMatthews & UsonI feOO^) and IlCriwattanawong et al 


1 201 ih . However, specially relevant is the comparison with proper 
Remagnitudes, whic h we highlight in Figure|3using coloured sym¬ 
bols for the works ofLHeraudeau_^^imien| l ll9^.l(^l de Paz et al.l 


l l2003l) . lKassin. de Jong, & Pogge 1 20061) and Amorln et alj 1 2007 ). 
The number of overlapping galaxies with Rc-band data is 
scarce, only 12 galaxies, but those are in good agreement with 
most previous works, specially with the more recent survey of 
iKassin. de Jong, & Poggd | |200^ . 

In the top pan el of Figured] we co mpare our isophotal radius 
with the results of Ijames et alj ( l2004l) . As the reference isopho¬ 
tal levels are different, there is a systematic offset in the isopho¬ 
tal radius in the sense that r esults from o ur work are systemati¬ 
cally smaller than those from Ijames et all but still there is a good 
correspondence between the two datasets. In the bottom panel of 
Figure [8] we show that the ellipticities are also well correlated, as 
expected, because at both reference isophotal levels the disk is the 
dominant component in the light of the galaxy, and has a simple 
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Figure 5. Difference between the surface brightness as function of the isophotal level for 12 galaxies profiles in common with de Jong and van der Kruit 
(1994) in t he Rc-band. The solid black (red ) lines represent the data from the OHP (SDSS), while the dashed lines are the profile errors. Multiple profiles of 
galaxies in ide Jong and van der KruiJ jl994) are shown in different lines. Dashed lines represent the errors in our profiles. 




“1^22.5 (mag) rso'’ (arcsec) (mag) 


Figure 6. Comparison of photometric properties derived for the 66 galaxies in common between the OHP and SDSS datasets. The parameter A represents the 
difference between the parameters in the two data sets, with its the mean value (< A >) and its standard deviation ((t(A)) shown in the bottom of each panel. 
Parameters are displayed as the following: (a) isophotal position angle; (b) isophotal ellipticity; (c) isophotal radius; (d) isophotal apparent magnitude; (e) 
effective radius; (f) total apparent magnitude. The horizontal and vertical axes display the results for the OHP and SDSS datasets respectively, and the dashed 
line represents the ideal equality between datasets. 
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i?23.5 (arcsec) - this work 



Figure 8. Comparison of our isophotal radius (above) and ellipticity (be¬ 
low) with those derived bv ljames^ HI <2oo 4) . Filled (hollow) points rep¬ 
resent data from OHP (SDSS), and the dashed lines represent the equality 
between measurements. 


geometry that does not vary drastically between these pass bands. 
Finally, in figure we compare the isophotal position angles and 
inclinations with the results from the Fla map analysis of Epinat et 
al. 2008, which demonstrate that our analysis produces results that 
are similar even to other tracers of the galaxy shape such as the gas. 


6 SCALING RELATIONS 

Scaling relations contain important information about the physical 
processes regarding galaxy formation and evolution, and impose 
impo rtant constraints to m odels that attempt to describe such ob¬ 
jects dCourteau et alj200% . In this section, we derive the most sig¬ 
nificant scaling relation involving luminosities, sizes and rotation 
curve velocity for each of the two most basic structural components 
of the galaxy, the bulge and the disk, and also for the whole galaxy. 
In section IfiTI we show how we correct the sizes and luminosities 
for the effects of distance, inclination and dust attenuation, and in 
section [ 6 ^ we show how we estimate the scaling relations. In sec¬ 
tion [63] we present the main results and in section I6A1 we explore 
the Tully-Fisher relation in greater detail. 




Figure 9. Comparison of our isophotal position angles (above) and inclina¬ 
tions (below) with those derived from the Ha velocity fields of Epinat et 
al. 2008. Filled (hollow) points represent data from OHP (SDSS), and the 
dashed line represents the equality between measurements. 


6.1 Correction for the effects of inclination and distances 

We use all galaxies in table|4|with a bulge component according to 
our classification in the decomposition. Disk apparent magnitudes 
(rridisk) are calculated using equation (O or by numerical integration 
in the case of broken profiles for the disks, while bulge luminosities 
(mbuigc) are calculated using equation { 8 |(. The absolute magnitudes 
of these components are then obtained using the equations 

AFdisk = m-disk — di — d 2 {l — cosi)'^^ — 51ogD — 25, (10) 

Mjuige = m-buige — hi — 62(1 — cosi)*'^ — 51ogD — 25, (11) 

where the internal extinction coefficients 61 = 0.60, 62 = 1.33, 
63 = 1.75, di = 0.15, d 2 = 1-09 and dz = 2.82 are obtained 
by linear interpolation from Table 1 of [Driver et alJ bOOSh for the 
Rc-band (A = 64 7nm), D is the distance in Mpc according to 
[Epinat et al.l ([2008[) . a nd the inclination i is taken from the gas ve¬ 
locity field analysis in [Epinat et al.[ l l2008h . if available, or from the 
isophotal analysis otherwise. Individual distance errors are rarely 
available, and we adopt a value of 25% for all objects. 

The total luminosity of each galaxy is estimated by its total 
magnitude according to the analysis of the curve of growth of the 
SB profiles (see section [Aa . In this case, we obtain the absolute 
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magnitude of the galaxies (Mb., total) from the apparent total mag¬ 
nitudes (ruR, total) using the expression 


Mr, total = mR, total “ Ai{R) - 5 logD - 25 , ( 12 ) 

where Ai(R) is the internal extinction correction given by 
iTullvetalJflg^ 

Ai{R) = log( 6 /a) X { 1.15 -h 1.88 (log 2 V„ax - 2 . 5 )} . ( 13 ) 

Here 6/a is again the minor-to-major axis ratio and Imax is the max¬ 
imum velocity of t he Hq rotation cur ve derived from the velocity 
field analysis from lEpinat et al.l l l2008h . 

We compare these luminosities with the physical sizes of each 
component, using the scale length of disks (h), or the inner disk 
length in the cases of broken disks, the effective radius of the bulge 
(re) and the effective radius of the galaxy (rso). We do not attempt 
to correct the sizes of the components for extinction, and we only 
rescale the sizes according to the dis tance. Finally, we use the ve¬ 
locity Imax as in lEpinat et al.l ( l2008ll as our dynamical tracer, ex¬ 
cluding galaxies for which the flat part of the rotation curve was 
not reached in the velocity field analysis of the Hq observations. 


6.2 Fitting method 

We assessed the statistical significance of pairs of luminosities, 
sizes and velocities using the Spearman’s rank correlation coeffi¬ 
cient r, which is a measurement of the strength of the correlation 
of two variables, and the associated p-value p, which indicates the 
probability of obtaining a result at least as extreme as the one ob¬ 
tained from a random distribution, both indicated in the boxes at 
the top of each panel. For 18 cases, we obtained correlations with 
p < 0.1%, for which we calculated scaling laws considering the 
direct and inverse cases. The relations, displayed in the form of 
dashed lines in Figure [Tol were calculated as in the following. We 
consider a linear relation in the form of 


yi = axi + P, (14) 

where each galaxy is represented by the index i, a is the slope of 
the relation and P is the zero point. We then performed a min¬ 
imization considering the measurement uncertainties in both vari- 
ables, considering al so an intrinsic scatter, eo, for the relation (see 
iTremaine et alj2003l . using the relation 
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1 
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iVi- P- axjf 
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(15) 


where N is the number of galaxies of the sample, v = N — 2is the 
number of degrees of freedom, and Eyi are the parameter uncer¬ 
tainties. The presence of the variables in both the numerator and in 
the denominator of relation (H makes the equation non-linear, and 
most common metho ds of minimizati on, such as the Levenberg- 
Marquardt algorithm jPress et al .Il992t) , may have problems to ob¬ 
tain convergence. To obt ain stable solutions, w e used the interac¬ 
tive method described bv ISedregal et alj l l2006h . which consists in 
solving equation GU for a fixed value of eq, and then update eq 
by multiplying for xt elevated to a power of 2/3, until obtaining 
xt = 1. This method failed in only two cases, as show in table 
[S] because either the dispersion is too low and x^ < 1 or if the 
dispersion was too high. The uncertainties for the coefficients were 
estimated using the bootstrapping method. 


6.3 Results 

Figure [Tol shows the resulting relations among luminosities, sizes 
and velocities for 62 galaxies selected in the previous section. We 
adopt different colouring for the panels above and below the diag¬ 
onal according to the strength of the bar and to the morphologi¬ 
cal classification respectively. However, we could not study these 
correlations in subsamples due to the low statistics after dividing 
the data in those classes, and our quantitative results are all related 
to the complete photometric sample. The summary of the scaling 
laws is shown in table[5] where we sort the relations by decreasing 
Spearman’s coefficients. 

All equations in table [5] can be used as a way of obtaining 
approximate physical parameters for one given measurement as 
well as for constraining models of galaxy formation at the cur¬ 
rent cosmic time. Out of the 21 pairs of parameters, we observe 
that only three combinations have relatively low correlation coeffi¬ 
cients. This indicates that most of the spiral galaxy properties are 
somehow linked. Although there are many possible properties that 
shape galaxies, such as different mass, angular momentum, and de¬ 
spite secular evolution effects, such as those which may form bars, 
there is still a great similarity among spirals which is still to be 
explained. Also, this large number of correlations restrict the inter¬ 
pretation of the correlations individually, and certainly a compre¬ 
hensible interpretation will be possible only wit h a more complete 
model of galaxy evolution (see lshen et al.ll20l^ . Nevertheless, we 
are going to briefly discuss a few of the scaling laws that have been 
observed here and previous ly in the literature, wit h the exception 
of the Tully-Fisher relation ^Tullv and Fishei[|T977h . which we ad¬ 
dress with greater detail in l6.4l 

The relation (q) between the sizes of bulges and disks 
was obtained previ ously by other authors dCourteau et alj[l9^ : 
lAguerri et alJuPOSil and may have impor tant clues for galaxy for¬ 
mation scenarios. ICourteau et ^ il l996h argue that disks should 
have been formed earlier than bulges and, therefore, the proper¬ 
ties of the bulges are linked to their host disks. Due to the rela¬ 
tively low Ser sic indices of the bulges, these are indeed likely to be 
pseudobulges dFisher and Drorvll200^ . which are formed by secu¬ 
lar evolution of the disks and, therefore, correlations among these 
parameters naturally arise in a scenario of secular evolution, disfa¬ 
voring scenarios of decoupled size relation such as bulges formed 
by mergers. We have found that the median value of Vo/h is 0.14 
considering all galaxie s, which is in agreement to the values in the 
literature (for instance. [Laurikainen et alj|2010[) . 

The luminosity of bulge is also of importance to understand its 
origin. Bulge luminosities and sizes are expected to correlate, as al¬ 
ready indicated in equation Lfauige oc Tg, and indeed there is a 
strong correlation as shown in equation (1). Moreover, the bulge lu¬ 
minosity is correlated to all other measured properties of the disks, 
to the whole galaxy and also to the rotation velocity. Therefore, 
the properties of the bulges we observe in late-type galaxies of the 
local universe are probably the result of secular evolution. Interest¬ 
ingly, the bulge lu minosity is also correlated with t he supermassive 
black hole masses ^Kormendy & Richstonall99^ . illustrating the 
important role of the bulges to understand the processes of galaxy 
formation yet to be fully understood. 

Another parameter that correlates strongly with almost all the 
others is the luminosity of the galaxy, as shown in equations (a), 
(b), (f), (g) and (j). The importan ce of the total luminosity, also 
observed bv ICourteau et alj <20071 ). may indicate that the baryonic 
portion of the galaxy has a pivotal role in the appearance of galax¬ 
ies: it is connected with the gravitational potential through the ve- 
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Morphology T-type Morphology 
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Figure 10. Rc -band correlations between luminosities and sizes of bulges (Mbuige^ log ^e), disks (Mdisk log h) and the galaxy luminosities (M^^total)’ 
sizes (rso) and maximum velocity of the rotation curve (Imax)- Panels above the diagonal display colours of galaxies according to the colour map in the 
bottom right of the figure, which divide galaxies according to the presence of bars following the morphological classification. Panels below the diagonal 
display colours according to the colourmap in the bottom left of the figure, classifying galaxies according to their numerical type in the de Vaucouleurs’s 
classification. The Spearman’s rank correlation coefficient (rg) and the p-value of the relation are presented in the box on top of each panel. For relations with 
significance greater than 3cr, we include a dashed black line to indicate the best fit linear regression, whose coefficients ai‘e presented in TableO The diagonal 
panels display the distributions of parameters as histograms. 
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Table 5. Scaling relations with statistical significance above 3a for the luminosities, sizes and velocity of the galaxies and its two basic subcomponents, bulges 
and disks, ranked by decreasing Spearman’s rank coefficients (r^). The first column indicates the identification of the relation. The second column indicates 
the variables involved in the relation, as well as the correlation coefficients and their p-values. The third column indicates the direct and inverse relations. The 
fourth column shows the intrinsic scatter of the relation. 


ID 

(1) 

Parameters 

(2) 

Relation 

(3) 

eo 

(4) 

(a) 

^R,totar^disk 
|r| = 0.96, p = 3 ■ 10-34 

Mdisk = (0.98 ± 0.03 )Mr -b (-0.6 ± 0.7) 

j^R,total = (1-02 ± 0.04)Mdi,j^ -b (0.6 ± 0.8) 

— 

(b) 

log k^max-Mj^ (Qjjj 
|r| = 0.88, p = 6 ■ 10-21 

Mr total = (—8.0 ± 0.7) log Vmax -b (—3 it 1) 
log imax = (—0.12 ± 0.01)M]^ -b (—0.2 it 0.2) 

0.37 + 0.05 

0.04 + 0.01 

(c) 

log ymax-Mjigj, 

|r| = 0.85, p = 2 ■ 10-1® 

Mdisk = (—7.1 it 0.8) log Vmax -t (—5 it 2) 
log Vmax = (-0.13 it 0.01)MHi,|, -t (-0.6 it 0.3) 

0.31 + 0.07 
0.04 + 0.01 

(d) 

log 7-50-log 

|r| = 0.85,p = 4 - 10-1® 

logh = (0.88 it 0.07) logrso -t (0.00 it 0.04) 
logrso = (0.93 it 0.06) log h-b (0.11 it 0.03) 

0.132 + 0.003 
0.135 + 0.004 

(e) 

log 7-50 

|r| = 0.83, p = 9 ■ 10-11 

Mdisk = (-4.7 it 0.3) logrso -t (-18.1 it 0.2) 
logrso = (-0.18 it 0.02)M^jij,j, -b (-3.2 it 0.3) 

0.61 + 0.03 
0.12 + 0.01 

(f) 

^R, total ‘-^bulge 
|r| = 0.82, p = 2 ■ 10-1® 

^bulge = (1-1 ± 0.1)MR mtal + (5 i 3) 
j^R,total = (0-65 it 0.06)Mbuige + (-9 it 1) 

1.11 + 0.02 

0.84 + 0.02 

(g) 

logrso-MRjQt^l 

|r| = 0.78, p = 9 ■ 10-14 

^R total = (-4-4 ± 0.5) logrso -t (-18.1 it 0.3) 
logrso = (-0.16 it 0.02 )Mr jqJ 31 -b (-2.7 it 0.3) 

0.88 + 0.02 

0.167 + 0.005 

(h) 

log Vmax-Mbulge 
|r| = 0.76, p = 7- 10-1® 

Mbuige = (—9 it 1) log Vmax + (2 it 2) 

log Vmax = (-0.08 it 0.01)Mbuige + (0.8 it 0.2) 

1.17 + 0.07 

0.107 + 0.007 

(i) 

log^-A7disk 

|r| = 0.72, p = 3 • 10-11 

Mdisk = (-4.3 it 0.4) log h + (-18.6 it 0.2) 
log h = (-0.16 it 0.02)Mdig|j + (-2.8 it 0.4) 

0.86 + 0.02 
0.165 + 0.006 

G) 

^R,total-l°g^ 

|r| = 0.72, p = 5 • 10-11 

log/i = (-0.15 it 0.02)MR mtal + (-2-5 i 0.3) 

MR total = (-4.3 it 0.6) log h + (-18.4 it 0.3) 

0.181 + 0.004 

0.98 + 0.01 

(k) 

'^disk'^bulge 
|r| = 0.69, p = 4 • 10-1° 

^bulge = (1-0 ± 0.1)Mdisk + (3.0 it 3.0) 

Mdisk = (0.55 it 0.09)Mbuige + (-11.0 it 2.0) 

1.42 + 0.03 

1.03 + 0.02 

(1) 

logre-Mbuige 
|r| = 0.67, p = 3 • 10-°9 

^bulge = (-4-1 ± 0.6) log re -t (-19.6 it 0.1) 
log re = (-0.14 it 0.02)Mbuige + (-2.9 it 0.3) 

1.41 + 0.02 

0.262 + 0.004 

(m) 

log Vmax-log rso 
|r| = 0.65, p = 1 • 10-°® 

log rso = (1.1 it 0.2) log Vmax + (—1.8 it 0.4) 
log Vmax = (0.5 it 0.1) logrso + (1.95 it 0.06) 

0.204 + 0.007 
0.134 + 0.005 

(n) 

log Vmax-log h 
|r| = 0.63, p = 3 • 10-°® 

log h = (1.0 it 0.1) log Vmax + (-1.7 it 0.3) 
log Vmax = (0.46 it 0.06) log h + (1.98 it 0.04) 

0.205 + 0.008 
0.141 + 0.004 

(o) 

log /i-Mbuige 
|r| = 0.60, p = 3 • 10-°1 

Mbuige = (-4.6 it 0.8) log h + (-16.3 it 0.4) 
log/i = (-0.09 it 0.01)Mbuige + (-1.1 ± 0.3) 

1.66 + 0.02 

0.231 + 0.002 

(P) 

logrso-Mbuige 
|r| = 0.53, p = 1 • 10-°® 

Mbuige = (-3.7 it 0.8) log rso+ (-16.6 it 0.5) 
logrso = (-0.08 ± 0.02)Mbuige + (-0.8 ± 0.4) 

1.83 + 0.02 

0.260 + 0.004 

(q) 

log h-log re 

|r| = 0.48, p = 9 • 10-°® 

log re = (0.8 ± 0.3) log + (-0.6 ± 0.2) 
log h= (0.3 ± 0.2) log re + (0.62 ± 0.03) 

0.41 + 0.01 

0.28 + 0.01 

(r) 

log 7-50 -log Te 

|r| = 0.41, p = 1 • 10-°® 

log re = (0.6 ± 0.2) logrso + (—0.5 ± 0.1) 
logrso = (0.3 ± 0.1) log re + (0.68 + 0.02) 

0.362 + 0.002 
0.275 + 0.003 
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locity, but has a more direct link with the size (or shape) of the 
galaxy. 

Other photometric relations well documented in the litera- 
ture include the relation (k) between bulge and disk luminosities 
dLaurikai nen et al. ll201(]|), and (f) which relates bulge and total lu- 
minosities dCarollo et al.ll2007h . The rotation velocity of the galax¬ 
ies is usually studied in comparison with integrated photometry, 
such as given in relation (b), the Tull y-Fisher relation , and t he size- 
velocity relation (m) also studied by ICourteau et al.l d2007h . How¬ 
ever, here we show that the rotation velocity also strongly corre¬ 
late with the luminosity and size of the disk component, as shown 
in relations (c) and (n), which is expected because the disk is re¬ 
sponsible for the majority of the light of the galaxy. Interestingly, 
however, the luminosity of the bulge also correlates with the ro¬ 
tation velocity, as shown in relation (h), indicating that the bulge 
properties have a dynamical link with the galaxy that hosts it. 

6.4 Tiilly-Fisher relation 

The T ully-Fisher relation (hereafter TF relation, ITuIIv and Fisher! 
Il977h is the most important scaling relation for disk galaxies, and 
it has been used for several purposes including distance determina¬ 
tions and, historically, as a way of measuring the Hubble constant. 
The TF relation relates the maximum velocity of the rotation curve, 
Vmax with the total magnitude of the galaxies in the form of a power 
law. The TF relati on has already been measur ed in the GHASP sam¬ 
ple previously by iTorres-Flores et al] d201lh in the near-infrared 
bands H and K, so here we add to those results the optical Rc-band. 
We adopt the following parametrization 

Ma = aA-log(^) -f/3A, (16) 

where M\ is the absolute total magnitude in the passband A, Hmax 
is the maximum velocity of the rotation curve, ax and Px are the 
slope and the zero point of the TF relation. Both the slope and the 
zero point of the Tully Fisher relation are of importance because 
they may be used either as constraints or as tests for models of 
galaxy formation and evolution. 

To produce a suitable sample for this specific relation, we se¬ 
lect the galaxies according to the following criteria. We remove 
galaxies with inclinations greater than 75° due to their high in¬ 
ternal extinction, and also galaxies with inclinations smaller than 
20 °because of their higher uncertainty in the determination of the 
rotation curve velocity. We also exclude galaxies with recession ve¬ 
locities lower than 3000 km s~^ due to possible peculiar velocities 
affecting the Hubble flow, except for the cases where more accu¬ 
rate distance indicators were used, such as Cepheids and red-giant 
branch distances. Finally, as we are only dealing with Ha velocity 
fields, we use the analysis of lEpinat et alj 120081) to exclude from 
the sample galaxies for which the maximum rotation velocity is 
not achieved according to their classification of the maps. We use 
our two absolute magnitude estimators, the asymptotic (Mr, total) 
and the isophotal (Mr,23.5) as the probe of the galaxy luminos¬ 
ity, resulting in samples with 80 and 72 galaxies respectively. Most 
galaxies of GHASP sample are not part of clusters of galaxies, so 
we consider that our TF relation is basically probing the field en¬ 
vironment, although the expect ed difference of the TF relation in 
differ ent environments is mild dPe Riicke et al.ll2007l : IMocz et al.l 

[20H . 

The TF relation in the Rc-band is shown in Figure fTTI To cal¬ 
culate the regression coefficients, we have used the minimiza- 


Table 6. Regression coefficients for the Rc-band Tully-Fisher relation using 
the total asymptotic magnitude and the total isophotal magnitude including 
80 and 72 galaxies respectively. 


Relation 

O-R 


£0 

TF 

Total asymptotic magnitude 
-8.1 ±0.5 -3.0 ±1.0 

0.28 ±0.07 

Inverse TF 

-8.4 ±0.7 -2.0 ±2.0 

0.28 ±0.08 

TF 

Total isophotal magnitude 
-8.9 ±0.6 -1.0 ±1.0 

0.33 ±0.06 

Inverse TF 

-9.3 ±0.7 

0.0 ±2.0 

0.34 ±0.08 


tion of equation 05} as descri bed in sectionl6.3l Also, the so-called 
inverse Tully-Fisher relation l lSchechteilll98(ll) coefficients are cal¬ 
culated as follows. We calculated the coefficients a' and /?' by in¬ 
terchanging the variables Xi yt in equation d, and then cal¬ 
culated the inverse TF relation using the relations a = 1/a' and 
P = —P'ja'. The summary of the results for the TF relation and 
the inverse TF relations is presented in Table|^ 

The TF relation is not just an important tool for measur¬ 
ing distances of galaxies, but it is also crucial to highlight pro¬ 
cesses of galaxy evolution, for instance, by comparing the TF re¬ 
lation of different morphological types. Spiral galaxies have a sin¬ 
gle TF relation, but lenticular and elliptical galaxies have TF rela- 
tions which run approximately parallel whe n compared to spirals 
llBedregal et al.ll20(^ : IPe Riicke et al.ll200^ . We observe that the 
TF relation of spirals in the Rcband is well defined for almost all 
galaxies, with only two exceptions that are worth discussing. NGC 
12276, marked in the figure [TT] with a yellow halo, does not seem 
to have any special feature to be offset from the TF relation, so one 
possibility to explain its position is that the distance to the galaxy 
is not accurate. We have used the value of 78.8 Mpc from Epinat et 
al. (2008) for consistency with the previous works, which is the ex¬ 
pected value ac cording to the systematic v elocity using the Hubble 
flow. However, [Pedreros & Mador3 ( Il98lh have estimated the dis¬ 
tance to this galaxy of 40 Mpc using the ring size, which implies in 
a difference of «1.5 magnitudes that is enough to bring the galaxy 
much closer to the TE relation. The galaxy with a red halo in fig¬ 
ure [TT] NGC 4256, has a peculiar morphology of a single arm and 
an asymmetric rotation curve, which may be the cause to move the 
galaxy off the TF relation defined by relatively more relaxed spi¬ 
rals. In this case, star formation may have been triggered recently 
as a response to a gravitational field, resulting in a relatively lumi¬ 
nous object compared to the TF relation. 

The slope and the zero point of the Rc-band TF relation in 
our w ork are in agreement wi th those previously derived the liter¬ 
ature. ITuIIv and Pierce! ( l200d) . for instance, determined the slope 
and zero point for a sample of 115 galaxies in four nearby clusters 
with velocities derived from HI line widths, and by not consider¬ 
ing errors in both variables nor the intrinsic scatter of the relation, 
they have found or = —7.65 and /3 r = —4.3, which is consis¬ 
tent w ith our results within 2 sigmas. On the other hand. lVerheiienI 
l l200lh has derived the TF relation for the Ursa Major cluster using 
the inverse relation without intrinsic scatter, and fixing the uncer¬ 
tainties in 0.05 mag for the magnitudes and 5% in the velocities. In 
this framework, they found slopes ranging from -7.1 to -9, and zero 
points ranging from -3.15 to 2.81 for their various samples, which 
is similar to our inverse TF relation parameters. 

The TF relation in the infrared pass bands is important be- 
cause these wavelengths a re reliable tracers of the stellar masses. 
iTorres-Flores et al] ( 1201 ih have used the GHASP sample to derive 
the TF relation in the H and K bands using 2MASS survey data 
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Figure 11. Tully-Fisher relation for the Rc-band for GHASP galaxies considering (a) the total asymptoti c magnitude , Mr. totals and (b) the isophotal total 
magnitude inside the isophote of 23.5 mag arcsec“^. The gray line in panel (a) indicates the results from ISakai et al.l l2Q0(Jl . In the top panel of each figure, 
the black dashed line indicates the direct TF relation results and the best fit while the red dashed line indicates the inverse TF r elation results. The 

bottom panels indicates the colour of the galaxies using the K band total magnitudes used for the TF relation in iTon'es-Flores et ii H2onh inside the isophote 
of 20 mag arcsec“^. The colour of each galaxy reflects the colourbar scale in the upper panels, separating objects according to their numerical morphology. 
The two objects highlighted with a halo are NGC 12276 (yellow) and NGC 4256 (red), which for different reasons are offset to the TF relation (see text for 
details). 


dSkrutskie et al.l[2005) as well as stellar and baryonic TF relations, 
and using a method similar to ours, they have obtained slopes of 
qh = —10.84±0.61 and qk = —11.07±0.63 and zero points of 
Ph = 1.97± 1.36 and Pk = 2.27± 1.39 for the H and K bands re¬ 
spectively. As expected, the R c-band slope is g reater than the slope 
in the near-infrared band (e.g. IVerheiienll20(^ . However, one im¬ 
portant feature observed in the infrared is a break in the TF relation 
for galaxies with log Vmax < 2.2, in the sense that galaxies below 
this velocity are under luminous related to the expected TF rela¬ 
tion for bright galaxies. This break in the TF relation is not noticed 
in the Rc-band. This difference in the shapes of the near-infrared 
and optical TF at the low-mass regime can be understood if one 
inspects the bottom panels of Figure[TT] which show the optical to 
near-infrared colours of the galaxies as a function of the maximum 
rotation velocities. These panels show a flat colour distribution ex¬ 
cept for the galaxies with Vjnax < 125 km s“^, indicating differ¬ 
ent mass-to-light ratios these galaxies. These low mass systems are 
bluer than more massive galaxies, indicating younger objects that 
may have been forming stars recently, and this effect fortuitously 
compensates for the difference in the stellar mass to light ratios, 
causing the differences in the shapes of the near-infrared and opti¬ 
cal TFs at low masses. 


7 SUMMARY AND CONCLUSION 

This study provided photometrically calibrated surface brightness 
profiles of GHASP galaxies, with 170 profiles in the Rc-band and 
108 in the SDSS bands u,g,r,i and z. From these data, we derived 


Rc-band integrated photometric parameters, presented in Table 
which are consistent with other works in the literature. All these 
results are public and will be available in digital format at the Fabry 
Perot repository in http: //cesam. lam. f r/fabryperot 

We perform multi component structural decompositions in the 
Rc-band, presented in Table [4] with the goal of separating the 
disk component from bulges, bars, lenses and nuclear sources, as 
a preparation to our forthcoming paper on the kinematic decompo¬ 
sition of GHASP velocity fields, which will be compared with the 
photometric work. 

Finally, we have applied new photometric data to observe 
bulges, disks and global scaling relations among luminosities, sizes 
and velocities in the Rc-band. We derived expressions for 18 scal¬ 
ing relations, which may be used to constrain models of galaxy 
formation and evolution. In particular, we studied the Tully Fisher 
relation using velocities derived solely from Ha maps from GHASP. 
We have obtained slopes and zero points that are consistent with 
previous findings in the literature in the Rc-band for cluster galax¬ 
ies, reinforcing the idea that the Tully Fisher relation is basically a 
relation between the total stellar content and the gravitational po¬ 
tential, which is barely affected by the environment and the pres¬ 
ence of photometric substructures. 
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